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Abstract 

Complex frequencies associated with quasinormal modes for large Reissner- 
Nordstrom Anti-de Sitter black holes have been computed. These frequencies have 
close relation to the black hole charge and do not linearly scale with the black hole 
temperature as in Schwarzschild Anti-de Sitter case. In terms of AdS/CFT corre- 
spondence, we found that the bigger the black hole charge is, the quicker for the 
approach to thermal equilibrium in the CFT. The properties of quasinormal modes 
for I > have also been studied. 
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Quasinomal modes of black holes have been an intriguing subject of discussions for the 
last few decades. It has become an evidence that the quasinormal ringing will dominate 
most processes involving perturbed black holes. This means that quasinormal modes 
will carry a unique fingerprint which would lead to the direct identification of the black 
hole existence. Detection of these quasinormal modes is expected to be realized through 
gravitational wave observation in the near future. In order to extract as much information 
as possible from gravitational wave signal, it is important that we understand exactly 
how the quasinormal modes behave for the parameters of black holes in different models. 
For black holes in asymptotically flat spacetime, especially spherical cases, they have 
been studied extensively, for a review see [1]. The study for the nonspherical black 
holes is developing [2,3]. Considering the case when the black hole is immersed in an 
expanding universe, the quasinormal modes of black holes in de Sitter space have also 
been investigated recently [4,5]. It was found that there are qualitative differences from 
the asymptotically flat case, in particular the scalar field decay is always exponential, 
rather than a power-law tail in asymptotically flat spacetime. This result has also been 
uncovered by a very recent study of the Schwarzschild AdS black hole model [6]. These 
observations support the earlier argument by Ching et al. [7] that usual inverse power-law 
tails as seen in asymptotically flat black hole spacetime, are not a general feature of wave 
propagation in curved spacetime. 

Motivated by the recent discovery of the AdS/CFT correspondence, the investigation 
of the quasinormal modes of AdS black holes becomes more appealing nowdays. The 
quasinormal frequencies of AdS black hole have direct interpretation in terms of the dual 
conformal field theory (CFT). In terms of the AdS/CFT correspondence [8-10], a large 
black hole corresponds to an approximately thermal state in the field theory, and the 
decay of the scalar field corresponds to the decay of perturbation of the state. After 
computing the scalar quasinormal modes of Schwarzschild AdS black holes in four, five 
and seven dimensions, Horowitz and Hubeny claimed that for large black hole both the real 
and the imaginary parts of the quasinormal frequences scale linearly with the black hole 
temperature. The timescale for approaching to the thermal equilibrium is determined by 
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the imaginary part of the lowest quasinormal frequency and is proportional to the inverse 
of the black hole temperature. However, for a small black hole the results are no longer 
the same as that for a large black hole. The quasinormal frequencies do not continue to 
scale with temperature. Some comments on small black holes were presented therein. 

The Schwarzschild AdS black hole studied in [6] is the simplest model in Anti-de Sitter 
space, which is determined by only two dimensionful parameters, the black hole event 
horizon r + and the AdS radius R relating to the cosmological constant by A = —3/R 2 . 
It is of interest to generalize the study of ref. [6] to a more general model, say Reissner- 
Nordstrom (RN) AdS black holes. Besides r + and R, the RN AdS black hole has another 
parameter, the charge Q. Thus it possesses richer physics to be explored. In this paper 
we are going to study RN AdS black hole in the hope of getting more understanding of 
how the quasinormal modes depend on this additional parameter and whether there is 
more information on AdS/CFT correspondence. 

The RN black hole solution of Einstein's equations in free space with a negative cos- 
mological constant A = —3/R 2 is given by 

ds 2 = -hdt 2 + h~ l dr 2 + r W, A = Q/rdt, (1) 

with 

r R z r r + r r z R z 
The asymptotic form of this spacetime is AdS. There is an outer horizon located at r = r+. 
The mass of the black hole is 

1 r 3 O 2 

The Hawking temperature is given by the expression 

r 2 R 2 



and the potential by 

= 

r + 



In the extreme case r + , Q satisfy the relation 

Q 2 3rl 

l-| + ^-0. (6) 

In the following discussions, we will concentrate our attention on the large black hole 
with r + ^> R. We will not consider the case of the small black hole here, partly because it 
is unstable, having negative specific heat [11] and partly because it is not of direct interest 
for the AdS/CFT correspondence [6]. 

Let us consider a massless scalar field $ in the RN AdS spacetime, obeying the wave 
equation 

□$ = (7) 

where □ = g a/3 V a V 'p is the d'Alembertian operator. If we decompose the scalar field 
according to 

$ = £-^r)*WM) (8) 

Ira T 

then each wave function ipi(r) satisfies the equation 

d 2 ipi d 2 ipi 



dt 2 dr* 2 



Vnl>i (9) 



where 



r z r dr 
JJLtH i r + + rl/R 2 + Q 2 /r + 2Q 2 2 

~~ a y r 2 + r 3 r 4 + R 2i 

dr 

and r* here is the tortoise coordinate defined by r* = J —. 

h 

The potential U has the same characteristic as that in Schwarzschild AdS black hole. 
It is positive and vanishes at the horizon, however it diverges at r = oo, which requires 
that $ vanishes at infinity. This is the boundary condition to be satisfied by the wave 
equation for the scalar field in AdS space. 

Quasinormal modes of AdS space are difined to be modes with only ingoing wave near 
the horizon. There only exists a discrete set of complex quasinormal frequencies. The 
quasinormal modes behave like e - laJ (*+ r *) near the horizon in RN AdS background. We 



thus introduce the ingoing Eddington coordinates by setting v — t + r*. The metric (1) 
can be rewritten as 

ds 2 = -hdv 2 + 2dvdr + r 2 {d0 2 + sin 2 #d0 2 ) (11) 

We separare the scalar field in a product form as 

$ = -i/;(r)Y(e,<l>)e- iwv (12) 
r 

The minimally-coupled scalar wave equation (7) may thereby be reduced to an ordinary, 
second order, linear differential equation with the radial terms yielding 

h(r)^Pp- + (h'(r) - 2iu)^P- - V(r)^(r) = 0, (13) 
ar A ar 

where the potential function is given by 
V(r) = ^ + ^±I> 

= i ( r± + A. + QL _ 2 Q1 + + <1+A. (14) 

ry ^"2 ^j^2^"2 ty ^ry^ ry^i J^j^ ^"2 

Note that by setting Q 2 = 0,R = 1, Eqs. (13,14) go back to the 4D Schwarzschild AdS 
black hole case addressed in [6]. In the following discussion we adopt R = 1. 

To find the complex values of uo such that (13) has a solution with ip finite at the 
horizon r = r+, and vanishing at infinity, we have to count on the numerical calculations. 
By using the numerical method suggested in [6] to compute the quasinormal modes, we 
will expand the solution in power series about the horizon and impose the boundary 
condition that the solution vanishes at infinity. Adopting the new variable x = 1/r, (13) 
can be reexpressed as 

+ v^ii^ + <^b** (a ° = (15) 

where the coefficient functions are 

s(x) = 16 

x — x + 

t(x) = 3r x 4 - 2x 3 - 4QV - 2itox 2 (17) 

u{x) = (x-x + )V(x) (18) 



and the parameter r = — - — + j~ — — —. As done in [61, we can expand s,t and u 

x% 

about the horizon x = x + in the form s(x) = J2 s n{x — x + ) n etc. The first terms are 
so = 2x\k, to = 2x 2 + (n — iu) and uq = 0, where k is the surface gravity and has the form 
k = (x + + 3/x + — Q 2 x\)/2. The solution of (15) can then be expressed as a power series 

oo 

1>(X) = J2 a n(x-X + ) n (19) 
n=0 
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Figure 1: Lines from the top to the bottom correspond to Q 2 = 10 7 , 10 6 , 10 5 etc. 

Substituting (19) into (15) and equating coefficients of (x — x + ) n for each n, we have 
the recursion relations for a n 

Y n—l 

a n = -— [K k ~ l)Sn-fc + ktn-k + Mn-fc]«fe (20) 
fe=0 

where 

P n = n(n — l)so + n ta — 2x 2 + n{nK — iu). (21) 

This relation has the same form as that of the Schwarzschild AdS case, but with an 
additional parameter Q herein. 
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r+ 

Figure 2: Lines from the top to the bottom correspond to Q 2 = 10 4 , 10 5 , 10 6 , 10 7 etc. 

In order to find the quasinormal modes for the AdS spacetime, we must select solutions 
which satisfy the boundary condition ip = as r — > oo(x — > 0). Thus we need to 
look for the zeros of Eq(19) at x = in complex uj planes. For a given l,x + ,Q the 
algorithm to find these frequencies follows the simple steps: (i) truncate (19) at a number 
iV of terms, construct the polynormial equation of uj using (20,21) so that (19) reduces 
to J2n=o a n{ UJ ){— x + ) n = 0, (ii) find roots of interest of this function, (iii) increase N 
until these roots become constant within the desired precision. Since the problem is 
to find numerical solutions of a polynomial equation, it becomes easy to use a built-in 
Mathematica function to locate zeros of J2n=o a n(^)(— x +) n directly. This procedure may 
reduce tedious trials and make the numerical calculations neat. 

We decompose the quasinormal frequencies into real and imaginary parts in the form 

uj = u r — iuji (22) 

which makes uJi positive for all quasinormal frequencies. For large black holes (r + 3> R), 
the relation of the values of the lowest quasinomal mode frequencies for I = and selected 
values of r + for different charge Q are exhibited in fig.l and fig. 2. The dots represent the 
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Figure 3: Lines from the top to the bottom correspond to Q 2 = 10 4 , 10 5 , 10 6 , 10 7 etc. 

lowest modes. In fig.l, lines from the top to the bottom correspond to Q 2 = 10 7 , 10 6 , 10 5 
etc. However the lines shown in fig.2 are corresponding to Q 2 = 10, 10 2 , ..10 7 from the 
top to the bottom, respectively. It is easy to see that with an additional parameter, the 
charge Q, neither the real nor the imaginary part of the frequency is a linear function of 
r + as found in Schwarzschild AdS case [6]. The bigger the charge Q is, the larger is the 
deviation from the linear relation we observe. 

The imaginary and real parts of the quasinormal frequencies relate to the damping 
time scale (t\ = l/c^) and oscillation time scale (r 2 = l/uj r ), respectively. From fig.l we 
learn that as Q increase, a>j increases as well, which corresponds to the decrease of the 
damping time scale. According to the AdS/CFT correspondence, this means that for big 
Q, it is quicker for the quasinormal ringing to settle down to the thermal equilibrium. 
Fig.2 tells us that the bigger charge Q leads to the smaller u r , which means that the 
frequency of the oscillation becomes small as Q increases. Therefore from fig.l and fig.2 
we can have a picture that if we perturbe a RN AdS black hole with high charge, the 
surrounding geometry will not "ring" as much and long as that of the black hole with 
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Figure 4: Lines from the top to the bottom correspond to Q 2 = 10 4 , 10 5 , 10 6 , 10 7 etc. 



small Q. It is easy for the perturbation on the highly charged AdS black hole background 
to return to thermal equilibrium. This is the new physics brought by the additional 
parameter Q in RN AdS black hole. 

^From (4), we learn that the temperature of the large charged AdS black hole does 
not scale linearly with the event horizon as that in Schwarzschild AdS black hole case. 
The behavior of the temperature is shown in fig. 3. The linear relation between T and 
r + is broken as Q increases. The relations between the real and the imaginary parts of 
the frequency and the temperature for large RN AdS black hole are displayed in fig.4,5. 
Again we see that in contrast to the results in [6], the relations are no longer linear when 
the additional parameter Q is taken into account. 

We have so far discussed only the lowest quasinormal modes with 1 = 0. Increasing 
I, we obtain the surprising effect of increasing the damping time scale (u>i decreases), 
and decreasing the oscillation time scale (u r increases) as in the case of Schwarzschild 
AdS black holes. As pointed out in [6], here we may also meet the problem of possible 
negative u>i as it continuously decreases with I. However from fig. 6, it looks that for the 
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Figure 5: Lines from the top to the bottom correspond to Q 2 = 10 7 , 10 6 , 10 5 etc. 

large black hole, the problem is not as serious as that shown in [6] for small Schwarzschild 
AdS case. Despite the similar behavior, once again the additional parameter Q led us to 
further new properties. As fig. 6 shows, different values of Q do not change the qualitative 
characteristic of decreasing Ui with I. Their decreasing rates are the same and do not 
depend on the value Q. However it is clear that the higher is the charge of AdS black 
hole, the later we will confront the tough question of cjj as I — > oo. The dependence of uj r 
on I for different values of Q are exhibited in fig. 7. Lines from upper to the bottom are 
Q 2 = 10 3 , 10 4 ...10 6 , respectively. As shown for the case of the imaginary part, the lines of 
different charges are parallel. 

In our numerical calculations, we found that we need for a large number N of terms in 
the partial sum to reduce the relative error in the computation of quasinormal frequencies 
as r + decreases and Q increases. However when the charge Q increases to nearly the ex- 
treme value satisfying (6), we cannot accurately determine the mode, no matter how large 
is the value of N that we adopt. The numerical convergence problem can be attributed 
to the method we adopted to compute the quasinormal modes. In the expansion of the 
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Figure 6: Dependence of Ui on I for large black hole with r + = 100. Lines from the top 
to the bottom on the left figure correspond to Q 2 = 10 5 , 10 4 and 10 3 . For the right figure, 
from the top to the bottom Q 2 = 10 7 , 10 6 , 10 5 , 10 4 . 



solution of the differential equation into power series, we have the convergence radius 



li m I ^ n I ^ I 3 + x 2 + Q 2 x A Sr 

■ Vn+i 9 - 3x+ + Ax\ - x% - Q 2 {hx% - x 5 ^ 1 



To expand the solution in a power series about the horizon, we need L > x + to ensure 



9 - 3x+ + 4x 2 , - x 



that the expansion is valid. In the range < Q 2 < % —, L increases from 

E\ r \ '' ' 

Lq = (3 + x 2 + ) /(9 — 3x + + Ax 2 + — x\ ) to infinity, where L > x + for a big black hole, which 
means that L > x + always hold in this range. However for Q 2 > ^-r— i — , L 



3 + x 2 

continuously decreases with increasing values of Q. Finally when Q 2 — > -r— - = 3r 4 +ri, 

x\ 

which is the exact extreme value, L — ► < x+. Therefore although the numerical 
approach is very efficient for determining the quasinormal modes for Schwarzschild AdS 
black hole and lowly charged AdS black hole, it breaks down for the nearly extreme RN 
AdS black hole case. 

It is important to notice that Eq(15) has the same form as the generalized spheroidal 
wave equations which also arises both in the quantum scattering theory of nonrelativistic 
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Figure 7: Dependence of oo r on I for large black hole r + = 100. Lines from the top to the 
bottom correspond to Q 2 = 10 3 , 10 4 , 10 5 , 10 6 etc. 

electrons from polar molecules and ions and in the theory of radiation process involving 
black hole in asymptotically flat spacetime. Similar convergence problem also appeared 
and hindered the earlier study of the quasinormal modes for extreme charged black hole 
and extreme rotating black hole [11,12]. 

Besides problems related to the method, is there some deep physics to account for 
the similar no convergence problem in determining quasinormal modes for extreme black 
holes both in asymptotically flat and AdS spaces? We know that the frequencies and 
damping times of the quasinomal modes are entirely fixed by the black hole, and are 
independent of the initial perturbation. It has been shown that there is a second order 
phase transition in the extreme limit of black holes [14-16]. This result has been further 
promoted in a recent study for the charged AdS black holes [11,17]. Because of the phase 
transition, the fluctuation of thermodynamic quantities become tremendous. It is hard 
to expect that we can obtain the fixed quansinormal frequencies to characteristize the 
thermalization timescale in the strong coupled CFT on the extreme RN AdS background. 
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We speculate that the black hole phase transition maybe the candidate physical reason 
behind this problem. 

In summary, we have computed the scalar quasinormal modes of large RN AdS black 
hole. These modes govern the late time decay of a minimally coupled scalar field, such as 
the dilaton. Compared to the Schwarzschild AdS black hole, we found that the additional 
parameter Q in RN AdS black hole has brought in some new properties. We observed 
that these modes no longer linearly scale with the black hole temperature. By AdS/CFT 
correspondence, we can interpretate the decay of the quasinormal modes to the time scale 
approaching thermal equilibrium in CFT. We learnt that perturbations are associated 
with the black hole charge. The larger the charge of the RN AdS black hole, the sooner 
it returns to the thermal equilibrium. The dependence of the quasinormal frequencies on 
nonzero angular momentum I has also been discussed. Due to the convergence problem, 
the method we adopted cannot be extended directly to study of the extreme black holes 
cases. Further refinement of the numerical method of solving the differential equation 
(13) is called for. Based upon some critical phenomena uncovered, we have given some 
speculation of the possible physical reason behind the problem for the extreme black holes. 
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